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We investigate vector contributions to the Lagrangian of A3— massive gravity in the decoupUng 
limit, the less explored sector of this theory. The main purpose is to understand the stability of 
. - - maximally symmetric vacuum solutions. Around self-accelerating configurations, vector degrees of 

freedom become strongly coupled since their kinetic terms vanish, so their dynamics is controlled by 
higher order interactions. Even in the decoupling limit, the vector Lagrangian contains an infinite 
number of terms. We develop a systematic method to covariantly determine the vector Lagrangian 
at each order in perturbations, fully manifesting the symmetries of the system. We show that, 
around self-accelerating solutions, the structure of higher order p-form Galileons arise, avoiding the 
emergence of a sixth BD ghost mode. However, a careful analysis shows that there are directions 
along which the Hamiltonian is unbounded from below. This instability can be interpreted as one of 
, the available fifth physical modes behaving as a ghost. Therefore, we conclude that self-accelerating 

configurations, in the decoupling limit of As-massive gravity, are generically unstable. 
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I. INTRODUCTION 

One of the deepest problems in theoretical cosmology is to provide a satisfactory explanation for the observed 
acceleration of our universe. It is intriguing to ask whether this phenomenon may be due to a non-standard dynamics 



> 

] of gravity at large distances. Einstein's General Relativity cannot lead to cosmological acceleration without the 
\^ . addition of an energy momentum tensor of the suitable form. The simplest choice is a positive cosmological constant 
■ term, whose small value is, however, unnatural from an effective quantum field theory perspective. Modifications 
, to General Relativity might achieve acceleration by exploiting self-interactions of gravitational degrees of freedom, 
I ■ without adding an energy momentum tensor. This phenomenon is dubbed cosmological self-acceleration. 
] An early explicit realization of this idea is the model of Dvali, Gabadadze, and Porrati (DGP) T]. In this set-up, 
T-H ■ a brane is embedded in a higher dimensional space-time with extra-dimensions of infinite volume. Four dimensional 
IL^ I gravity is recovered at small scales thanks to an Einstein-Hilbert term localised on the brane. Above a cross-over 
scale, which depends on a combination of the four and higher dimensional Planck masses, gravity becomes higher 
dimensional, and this leads to cosmological self-acceleration. This scenario is appealing for its simplicity, and for its 
H . clear geometrical interpretation. However, the dynamics of cosmological fluctuations around self-accelerating solutions 
- - ■ is plagued by ghost instabilities; some of the dynamical modes have kinetic terms with the wrong sign [2|. A way out to 
this conclusion can be achieved by modifying and generalizing this set-up to develop scenarios with the right features 
to avoid the emergence of a ghost. A possibility is to adopt the Galileon proposal Q. An interesting model, so-called 
A3 massive gravity has been developed based on this suggestion, and admits self-accelerating configurations, 
with the scale of acceleration set by the graviton mass. This is a generalization of Fierz-Pauli massive gravity with 
the addition of higher order graviton self-interactions that are free of the Boulware-Deser 'sixth' ghost mode, which 
plagues the original Fierz-Pauli theory ^5;]. Around a flat Minkowski background, the theory propagates flve healthy 
degrees of freedom: two helicity 2, two helicity 1, and a scalar mode. However, one (or more) of these five modes may 
become a ghost around solutions that are different from Minkowski space. 

In the so-called decoupling limit Q, the structure of self-accelerating configurations is relatively easy to study, and 
the theory exhibits the explicit Galileon symmetry in the scalar sector, resulting in an attractive model from the point 
of view of stability under classical and quantum corrections (see 0,3 for reviews). If instabilities are found in the 
regime well described by the decoupling limit, the theory is certainly sick. On self-accelerating cosmological solutions, 
it has been shown that vector degrees of freedom become strongly coupled, being characterized by vanishing kinetic 
terms Our purpose is to analyze in detail the Lagrangian controlling these vector modes, the less explored sector 

of the theory so far, and to understand whether they can lead to instabilities around cosmological configurations. Even 
in the decoupling limit, the Lagrangian for vector degrees of freedom contains an infinite number of terms. Extending 
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the approach introduced in [lO|, we are able to develop a systematic method that allows to determine in a covariant 
way the vector Lagrangian at each order in perturbations, fully manifesting the symmetries of the system. In order 
to study the role of vectors in the dynamics, one can proceed in different ways. In we showed that the inclusion 
of a background vector, while maintaining the qualitative properties of the background solutions, removes the strong 
coupling of vector perturbations, but forces one of the five degrees of freedom to become a ghost thus ruling out self- 
accelerating solutions with a background vector switched on. More recently, a qualitatively similar result was obtained 
away from the decoupling limit, by turning on background anisotropy and considering Bianchi-type solutions jTl|. In 
this work, instead, we focus on solutions with no background vector field and analyze the dynamics of perturbations. 
Scalar-vector interactions, appearing in the Lagrangian at order higher than two in perturbations, are the key points 



for our analysis. At third order in perturbations, these couplings form the so-called p-form combination 12| (see also 
[13}). and lead to equations of motion with at most two time derivatives. The same structure continues at fourth 
order in perturbations, after a suitable field redefinition has been performed. Besides showing this explicitly at fourth 
order, we provide strong arguments indicating that the same pattern continues at higher orders in perturbations; this 
confirms from a new perspective that the theory does not contain a sixth ghost BD mode around self-accelerating 
solutions. On the other hand, a Hamiltonian analysis of the resulting system including higher order perturbations 
reveals that the Hamiltonian is unbounded from below, signaling an instability which in turn can be understood as a 
ghost within the available five propagating modes. Therefore, our analysis shows, by a direct study of the Lagrangian 
structure at higher order in perturbations, that self-accelerating solutions in A3 massive gravity are indeed unstable. 

Our analysis of the vector Lagrangian in the decoupling limit is sufficiently simple that allows to physically under- 
stand the origin of the instabilities around self-accelerating configurations. At the same time makes manifest that 
the vector Lagrangian in this theory organizes in a very non-trivial way at higher order in perturbations, avoiding 
the appearance of a ghostly Boulware-Deser 'sixth mode' by forming combinations corresponding to p-form Galileon 
combinations. This allows to appreciate from a new perspective how a theory (in this particular case A3 massive 
gravity) can subtly avoid Ostrogradski-type instabilities [3, [l^ , though it leads to a different type of instabilities 
around self-accelerating backgrounds due to the lack of the vector kinetic term. 

This work is organized as follows. After starting with a brief review of A3 massive gravity in decoupling limit 
in section [III we consider in section IIIII the theory in the minimal set-up, — a4 = (see the next section HI] for 
a definition of these parameters), in the presence of a bare cosmological constant. We are able to classify all the 
maximally symmetric solutions of this theory, as a function of the cosmological constant. We show that no healthy 
maximally symmetric solutions de Sitter or anti-de Sitter solutions exist in this case. In section IIVI we generalize 
our discussion to the full theory with arbitrary and 04, reaching the same conclusions. Section |V] is devoted to a 
summary of our results, followed by three technical appendixes. 



II. PRELIMINARIES: A3 MASSIVE GRAVITY IN A DECOUPLING LIMIT 

The A3 theory of massive gravity Q is designed to remove the Boulware-Deser ghost by an uplift of a decoupling- 
limit construction of massive gravity, which does not propagates a sixth scalar degree of freedom. The resulting 
Lagrangian is then the following 

C = ^V^iR-ll{g)-12A,,), (1) 

where R is the Ricci scalar, and U{g) is a potential term for the graviton that will be explicitly shown later. We add 
a pure cosmological constant Acc to the original formulation but avoid any additional energy momentum tensor, 
since we focus on maximally symmetric vacuum configurations. These solutions are the simplest examples of cosmo- 
logical solutions in the theory which, nevertheless, exhibit a rich behavior with respect to cosmological fluctuations. 
Cosmological solutions in this theory, and in the related ghost-free bigravity models, have been extensively studied 
[qI. [I6I - 12II] . In order to write down the potential U{g), it useful to express the graviton in a diffeomorphism invariant 
way, by means of the Stiickelberg trick; hence writing the metric as 

5^11' = Vp^u + -77^ = Hf,„ + riap df,(l)"d^<j)'^ , (2) 
IVlp 

where the four Stiickelberg fields 0" transform as scalars under Lorentz transformations of the physical metric g^^, 
while transforming as four vectors from the fiducial metric rjap perspective. Greek indeces starting from /i indicate 
quantities evaluated in physical space, while greek indeces starting from a indicate quantities evaluated in the fiducial 
space. Notice that we expand the metric g^^ about a flat background 7]^ij in physical space. Using Helmoltz theorem 
in the fiducial space, the helicity one, Ap, and the helicity zero, tt, modes of the Stiickelberg fields 0" can be extracted 
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by the expansion 

A3 A3 

where 



(3) 



A3 = {mHlpif\ (4) 

corresponds to the cut-off scale for this theory. The powers of m and A3 in eq. ([3]) are included in such a way to 
provide, after a suitable diagonalization, canonical kinetic terms for the tensor, vector and scalar degrees of freedom 
around Minkowski space-time. In the remaining part of this section, we set the vector fluctuations to zero, but will 
return to study them in detail in the following sections. 

We are interested in studying the so-called decoupling limit of A3 massive gravity, which is defined as 

m->0, Mpi ^ 00 , A3 = fixed. (5) 

Moreover, in this limit, we demand that the cosmological constant Acc survives such that 

A/|,Acc ^ A^Ao , Ao = fixed , (6) 

with Aq of dimensions of mass to the fourth. The resulting theory describes well the physics up to energies of order 
A3, where new higher dimensional operators arise (suppressed by new scales larger than A3). In this work we focus 
on scales smaller or equal to A3. When convenient, we use the abbreviations 

n^,. = dfj,d^TT and 11 = 9^9^7r . 

In the decoupling limit, it is equivalent to raise and lower indices with the physical or the fiducial metric; both 
reduce to flat space since the effect of can be neglected in the limit of infinite Planck mass. In order to render 
the formulae less cumbersome, we set the scale A3 = 1; if desired, it is straightforward to re-instate it by means of 
dimensional arguments. Then the tensor iJ^^ appearing in eq. ([2|) reads (neglecting for the time being the effect of 
the vectors) 

iJ^. = + 2n^, - T^^^U^^n^p. (7) 
IVIp 

It is convenient to define a quantity 



IC-, = 5;- - Hi:, (8) 
which has the property that /Cuz^L =o ~ ^uiy- By means of this quantity, we consider the following potential [3] 

U = -w? [U2 + a3U3 + ai Ua] , (9) 

with 

U2 = (tr/C)2 - tr(/C2), 

Uz = (tr/C)3-3(tr/C)(tr/C2) + 2tr/C^ 

Ui = (tr/C)''-6(tr/C)2(tr/C2)+8(tr/C)(tr/C3) + 3(tr/C2)2-6tr/C\ 

where m has dimension of a mass, while and 04 are dimensionless parameters. By construction, this potential 
leads to the following ghost free Lagrangian for tensor and scalar modes in a decoupling limit j3] 

U,., . = -Ih^'^^f.Kp -3Xoh + h^^ (xW -I- + ^(3)j ^ (10) 

where E^j^ is the operator acting on Zap as 

£^!^Zap - (az^, - d^daZ^ - a^d^z^ + a^a^z^ - ry^.nz^ + v^.MpZ^'^) . (ii) 
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The expressions for the are given by 



(1 + Sag) 
1 



nn^. + -77^. (n^ - n;;n[;) 



(as + 3q;4) 



6n?.„ - en nl^ + sn^, (n^ - m n;;) 



(12) 



This Lagrangian is at most quadratic in h^i, and the scalar tt couples to /i^i, through higher derivative interactions, 
which lead to equations of motion with at most two time derivatives. Notice that the Lagrangian (|10p contains a linear 
term proportional to Aq h, which indicates that Minkowski space, with cj)" = cc", is not a solution of the background 
field equations. In order to remove this linear term, one can set the cosmological constant to zero (Aq = 0), or consider 
more general background configurations, as we will do in the next section. 



III. MAXIMALLY SYMMETRIC CONFIGURATIONS WITH og = 04 = 

In this section, we consider a generalization of the maximally symmetric configurations found in including a 
cosmological constant, and discuss the dynamics of fluctuations around them. We focus here on the case = = 0, 
postponing the analysis of the general case to section Hvl The novelty of our discussion is the detailed study of the 
vector fluctuations around cosmological solutions. These cosmological backgrounds are solutions in the decoupling 
limit, so one might wonder whether there exist maximally symmetric configurations of the full theory which reduce 
to the ones we study here. The answer is afiirmative; we explicitly showed this fact in [10] for the solutions found 
in [1^, [13], and by the same methods, it should be straightforward to investigate other configurations found in the 
literature as [l^ (see also [T9|). 



A. Properties of tensor and scalar fluctuations 

An important property of the decoupling limit ([5]) is that it does not require the scalar tt to be small. The Lagrangian 
given in (jlO[) is correct no matter what size the scalar tt is. In other words, the decoupling limit Lagrangian is not 
a perturbative expansion in field fluctuations, but a perturbative expansion in the small scales m and 1/Mpi, while 
keeping A3 flxed. This is different from standard perturbation theory, in which instead the quantities involved are small 
and only contributions up to a given order in perturbations are included. We can exploit this fact to straightforwardly 
determine maximally symmetric solutions in this set-up. Let us parameterize the h^j,u and tt degrees of freedom as 

hf^i, = — ^x'^ 111^1^ + hf_ii, , (13) 
^ = (i_£El.^2^^, (14) 

where a;^ = x^x^rj^i,, and cq is a constant parameter. is a quantity with dimension of mass to the cube (this is 
the quantity dubbed MpiH^ in [9], with of dimension mass squared; in their notation, the scale of is set by the 
graviton mass rri^). A non-trivial proflle for the vector could be included, but we have already studied the subject in 
detail in [10] (at least for spherically symmetric set-ups). By a suitable choice of cq and H, the previous expressions 
solve the equations of motion for the lagrangian ([TU]). so that the background metric and Stuckelberg fields would 
then be given by 

9^ = (1 - ^x^) V,. , (15) 
Co . (16) 



and h^i, and tt will be the quantities playing the role of our helicity two and helicity zero 'small' fluctuations around 
this background. 

The metric corresponds to a maximally symmetric space: (A)dS for (negative) positive H^. In the decoupling 
limit, the metric of a maximally symmetric space with vanishing spatial curvature can be reduced to the form (jlSp 
by appropriate change of coordinates; a background profile for the Stiickelberg fields (fT6|) allows us to solve the 
corresponding equations of motion as we will to discuss in what follows. 
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Plugging the parameterization (|13l) into the Lagrangian for tensor and scalar fluctuations, eq. ^TU\\ with 013 = a4 
and using the following useful relations, which are valid up to total derivatives when included in the action 



one finds 



yy^'' [77^, Dtt - 8^3, tt] = 247r , 



1 f , 
2 

-12i?2(3-2co)7r + ft^''X^^ 



- -T^h^^'^t^Kp - 3 [Ao - - (1 - Co) (2 - Co)] h 

6H^ ttDtt . 



(17) 
(18) 
(19) 



(20) 



The expression for X^^ is given by 



Xfii, — 



(3 - 2co)77M.n - (3 - 2co)n^. - + n^riA. + -r;^. - n^^K 



We can perform a field redefinition that decouples helicity 2 from helicity fields 

hf_,u h^y + (3 - 2cQ)rif_,yTT - . 



(21) 



(22) 



Then the kinetic terms for tensor and scalar are diagonalized, resulting up to total derivatives and terms independent 
of fields 

U,^,. = -^Vf^""''/ia^-3 [Ao-i/2_(i_co)(2-co)] (/i + 4(3 - 2co)7r + fflTT 



2 

-12H^{3 - 2co)7r 



^(3-2co)' 



ttDtt - ^ (3 - 2co) {dn fOn 



(23) 



Let us emphasize a crucial feature of the previous Lagrangian: the decoupling limit and diagonalization procedures 
provide a contribution to the kinetic term of the scalar fluctuations, which includes a term proportional to the 
curvature of space-time. This term is inherited from a third order interaction between tensor and scalar, which 
after shifting the metric like in eq. (|13p . is obtained by applying the relation (jl9p . Perturbing around curved space, 
the scalar component of the massive graviton acquires a kinetic term depending on space-time curvature. This was 
already noticed in Q , in the context of the Fierz-Pauli theory. 

We then observe that the Lagrangian (|23)) contains two tadpole terms that depend linearly on the fields. Therefore, 
in order to ensure that our configurations (|15p - (|16p are solutions of background field equations, we have to remove 
these tadpoles by setting their coefficients to zero. The conditions that we find are 



(1 - Co) (2 - Co) 
i/2(3~2co) 



(24) 
(25) 



The last of these conditions gives two branches of solutions. We will discuss these background solutions and their 
stability under cosmological perturbations in Section IIII CI after analysing the dynamics of vector fiuctuations in the 
next subsection. 



B. Including vectors fluctuations 

Here we discuss in detail the Lagrangian for vector fluctuations, the less explored sector of this theory. We will 
follow and develop further the procedure originally elaborated in [lOj , that we briefly review now. In Section UTJ we 
defined the tensor Hf^^ as (recall that we set A3 = 1) 

H^, = g^, - (5^ - - m d^A") (J^ - - m d,A^) . (26) 

From now on, we can neglect the helicity two fiuctuation hfj_^, since in the decoupling limit, h^i, does not couple to 
the vectors, hence all possible interactions and self-interactions are those studied in the the previous section. Recall 
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that we raise and lower indexes with 77^^. It is convenient to adopt a matrix notation for tensors. If B^^, and C^jy are 
two tensors, we write 



B = B,''; BC = B/C/ ; BC^ = B.^C^ ; tri? = B, 



Let us introduce the tensor 



Then one can write 



with {dA = d^.A") 



(27) 



- = A" P^' + P/ d'^Ap = dAP + PdA^, 



L 



2p 



dnA°'d''Aa = dAdA'^ 



(28) 



We will need to take the square root of M up to second order in an expansion in m, so it is convenient to write it as 

(29) 



M = P'^ -mLi+ L2, 



= p2 [l-mQi+m^Q2]' 



with 



h - Ip'^L2 - Ip-^LiP-^Li . 



Then 



VaI = P (1 -mQi +TO2g2) + mD + m'^E + 0{m^) 

with D and E two matrices satisfying the following equations (see [l^l for details) 

{D,P} = P[Qi,P], 

= [Q2, P] + [D, Qi] + {P, E} + D^, 

where [, ] and {, } are the commutator and anticommutator respectively. Taking traces, one finds 

trP) = , tiE = -itr {P-^D'^) . 



(30) 



(31) 



(32) 



(33) 
(34) 



(35) 



The previous formulae are correct up to second order in an m expansion, and this is sufficient for the decoupling limit, 
since higfher order terms vanish. 

The vector-scalar Lagrangian that can be condensed in the following expression 



= 7 



6 trVM + tiM - ( trVM] - 12 



(36) 



Plugging the previous formulae for M and and focussing on the part inside the parenthesis that is proportional 

to (term with m*' and are total derivatives or field- independent, and those with higher powers in m vanish in 
decoupling limit 771 —> 0), one finds [l3| 



Ca, = - { [trLa - (trPQi)^] + 2 (3 - trP) tr (PO2) + (trP - 3) tr (p-^D^) } 



(37) 



We are now interested in the dynamics of fluctuations around the maximally symmetric configurations discussed in 
the previous section. In order to achieve this, we perform the shift ((Ti)) to the scalar background configuration, that 
amounts to express P as 



p = (co 1 - n) 



(38) 
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Because the above tensor P is obtained by taking the square root of P^, we are forced to set cq > 0. We decide 
to calculate terms up to quartic in perturbations (in the hat quantities). The solution for the matrix equation p3p 
determining D starts quadratic in perturbations, and it reads at this order 



D ^ — 



n, {dA + dA^ 



(39) 



Calculating the remaining contributions and assembling them together, one finds, up to quartic order in perturbations, 
the following action coupling vectors to scalars 



— ( (3 - 2co ) trF^ + trn trF^ + 
^ 'tr(fl2F2) +tr(nFnF) 



(3 - 4co) 



4 Co 

(3 - 4co) 



2 eg 



Co 

1 

Co 



— tr [flF^ 
trntr(ni^2)| 



(40) 
(41) 



This Lagrangian exhibits both a Galileon symmetry in the scalar sector, and an abelian gauge symmetry in the 
vector sector. The simple and symmetric form for this scalar-vector Lagrangian will allow us to study the stability 
of cosmological vacuum solutions in a straightforward manner. Notice that the previous terms are only the firsts of 
an infinite series that characterizes the vector Lagrangian in the decoupling limit [loj . Our method allows one to 
systematically compute each order in perturbations, in a covariant way that makes manifest the symmetries of the 
system. These initial terms are, in any case, sufficient for our purposes (see however Appendix |^ for a resummation 
of all the contributions in a special case). 



C. Mciximally symmetric configurations and their stability 



Armed with the previous results, let us now discuss the properties of the configurations that solve the two tadpole 
equations ([M)) and (P5l) . and the the dynamics of fiuctuations around these solutions. As we anticipated, eq. ([25)1 
gives two branches of solutions (from now on, for simplicity, we will remove the hats from the scalar fluctuations): 



Branch I: Minkowski space. Choosing H — to solve eq. p5|) we obtain Minkowski space (the resulting configu- 
ration corresponds to the screening solution discussed in [§]). In this case, the second tadpole eq. gives 



Co 



3, 1 r- 



4Ao) 



(42) 



which in order to have a well defined square root, forces Aq > —1/4. Let us consider now the Lagrangian for vector 
fiuctuations, eq. (|4ip . that we rewrite keeping terms up to third order in fluctuations and expanding the traces: 



_ (3 - 2co) 

4 CO '^'"'^ 4 CO 



(3 - 4co 

Co 



(43) 



We notice that the kinetic term for the vector fluctuations is healthy when co lies in the interval < cq < 3/2. 
This implies that only the choice with minus sign in eq. (I42p gives a healthy kinetic term for the vector, and the 
cosmological constant Ao in this branch should be contained in the interval — l/4<Ao<2in order to avoid ghosts 
and have a well defined solution. The third order terms in the Lagrangian P5)) can be removed by the following field 
redefinition of the vector field 



F„ 



Co 



(3 - 2co) 



2(2co - 3) 



2co(2co-3)''^"^''^^^ 



(44) 



This redefinition, valid as long as co lies in the interval < co < 3/2, renders the vector Lagrangian canonical up 
to third order in perturbations. Similar redefinitions should allow to remove terms beyond the cubic order in the 
scalar- vector Lagrangian. 

Branch II: (A)dS space. This amounts to choose co — 3/2 to satisfy eq. Then, the second tadpole, eq. (|24p . 

implies the condition ~ ^{1 + 4Ao). By inspecting the kinetic term for the scalar in eq (1^^ . one then finds that, 
in order to avoid a ghost in the scalar sector, the parameter must be negative (corresponding to AdS space). 
However, in this case cq = 3/2 the vector Lagrangian is strongly coupled, since the vector kinetic term vanishes. The 
vector Lagrangian (|4T|) reads in this case 



^ [3 (trntri^^ - 2trnF2) - 2 (trtfF^ 



trntrBi^^ + trni^ni^)] 



(45) 
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As mentioned above, the kinetic terms for the vectors vanish; however, vectors become dynamical by couphng them 
with the scalar at third or higher order in fluctuations (this was already pointed out in Nevertheless, one 

should worry about higher derivatives in the equations of motion, since the previous Lagrangian contains contributions 
with two time derivatives in the scalar field tt. For systems coupling scalars with vectors, it is possible to find the 
combination that ensures that the equations of motion do not contain at all terms containing more than two time 
derivatives. It is a generalization of Galileon combinations which was explored in [13] and dubbed p-form Galileons. 
Up to fourth order in perturbations, the correct combination (without including higher derivatives in A^) is 



^p-form ^ [trntrF^ - 2trnF2] + bo jtri^^ ^5.^2 „ (trH) 



- 4trntrni^^ 



2trni^nF 



} (46) 



where ao and 6p are arbitrary coefficients. The above third order action with the aforementioned properties was 
presented in [12], while the fourth order one is as far as we know new. Comparing (|45|) with (j46l) we notice that 
while the third order action has the correct structure to avoid higher order time derivatives in the equation of motion, 
the fourth order Lagrangian does not seem to satisfy this requirement. However, a suitable field redefinition allows 
to recast eq. (j45p into a healthy form: we will discuss this technical point in Appendix \^ also providing strong 
arguments that indicate that the same pattern continues at higher order in perturbations. 



On the other hand, although our scalar- vector Lagrangian (|45|) does not lead to a propagation of a sixth ghost mode, 
it does generally lead to a ghost-like instability around self-accelerating configurations, in which the ghost is one of 
the available vector modes. We have already shown in 10] that, when turning on a non-trivial profile for the vector 
solving the equations of motion, the corresponding Lagrangian for perturbations around the resulting configuration 
acquires kinetic terms for the vector with the wrong sign. For completeness, we include a proof of this fact in Appendix 
IbI based on the Lagrangian (j45p expanded up to third order in perturbations. In the remaining of this section, we 
instead directly point out the instability by analyzing the hamiltonian associated with the Lagrangian obtained by 
combining the third order lagrangian contained in (1451) with the scalar kinetic term: 



C"'''''^ = -SiJ^TrDTT - - [D-kF^^F"" - 2d^,^TTF>'PF''p] 



(47) 



We choose, for simplicity, the gauge Aq = 0, diAi — 0; by doing a standard 3 + 1— decomposition, the previous 
Lagrangian reads 



C 



1 r 
3 



2Tr AiAAi + An Af - ij A,Aj 



(48) 



where the dots represent the terms without time derivatives, which we do not include since they do not play a role in 
the present discussion. The conjugate momenta to tt and Ai are 



-6H^ 7T 



1 



■A; A A 



7T A Aj + AtT Ai - TT, ij Aj 



In order to analyze the associated Hamiltonian, it is convenient to introduce the matrix 

TT 



Kij = AirSij 



(49) 
(50) 

(51) 



written in terms of spatial derivatives of the scalar fluctuation; our analysis depends on the properties of this object. 
First, let us suppose that the scalars fiuctuations are such that this matrix vanishes: Kij = 0. Then, we can easily 
invert the relations that define the conjugate momenta, and obtain 



3A,Ha. 
2 A, A A,; 



A, AA, 

which translates into the following Hamiltonian 

Hi 



n = - 



1 



12i?2 12^2 



H, 



H. 



9g2 a,Ua, 
A, AA, 



A,; AA„ 



3 A,n^. 

2 " A, AA, 



27 



A,Ha. 
A, AA,; 



(52) 
(53) 

(54) 
(55) 
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where the dots represent terms without momentum variables. The previous hamiltonian is Unear in Hj^; hence it 
is unbounded from below. Notice that this argument holds even in the limit in which vanishes. In conclusion, 
perturbations of the background self-accelerating solution, along the direction of scalar fluctuations such that = 0, 
admit unstable directions along which the system falls towards regions where the energy is unbounded from below. 

Similar conclusions hold for more generic Kij. Let us, for example, consider a Hij that is non- vanishing, and 
invertible. Then, after straightforward manipulations, one can show that the Hamiltonian can be written as 



In = (n.-AA.«-n,,)%n,.«-n., + ..., (56) 

where, again, the dots represent terms without momentum variables. It is not difficult to see that there are many 
unstable directions associated with this Hamiltonian. For example, make a choice for the vector /S.Ai, so that the 
scalar combination C = AA^K^-^AAj is non-vanishing and has a certain sign; for definiteness, the magnitude of 
AAi is chosen such that the denominator of the first term has the same sign of C. Accordingly, choose the vector 
H^. such that H^.k~-^H^^ has the same sign of C (for example, choose it in the same direction of the AA^). Then, 
choosing suitably the magnitude of H^r, it is possible to make one of the two terms in the previous hamiltonian 
arbitrarily negative ~ hence the Hamiltonian is unbounded from below. Other cases, such as the case in which is 
non-vanishing but not invertible, can be treated in a similar way. 

In conclusion, one generically expects that instabilities arise. There are many directions in the moduli space of 
fluctuations along which the energy is unbounded from below, and towards which the system can be driven into 
dangerous regions. 

To summarize, for each value of Aq it is possible to determine a unique configuration that belongs to the first 
or second branch (see figure [1]). Branch I are the screening solutions discussed in [9j, where the allowed range for 
the cosmological constant Aq is small, from —1/4 to 2 in order to avoid ghosts. The value Ao — —1/4 is special, 
being common to both branches. In this case, scalar and vector degrees of freedom are strongly coupled. This is 
the case studied in (22j in the context of Fierz-Pauli massive gravity, where new symmetries were observed in linear 
perturbations. It is not clear whether these symmetries survive to higher order in perturbation theory. However, from 
Lagrangian (|41l) one can see that in this common point, the cubic terms do not vanish, suggesting the disapperance of 
the symmetry beyond linear order in perturbations in A3 theory: in fact, as we have found above, these higher order 
contributions make the case H — unstable. For Aq < —1/4, the cosmological constant curves the space leading 
to an AdS configuration that is unstable since it is characterized by a ghost in the vector or scalar sectors. As we 
anticipated, the analysis of vector fluctuations is essential for characterizing completely the two branches and exhibit 
the instabilities. In the next section, we repeat the analysis using configurations with non-vanishing as and a^. 
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FIG. 1: Branches I (solid vertical line with H = Q) and II (solid horizontal line) solutions for as = 04 = 0. Realty configurations 
are shown in blue. The branch II of self-accelerating solutions has no kinetic term for the vector field, and there are directions 
in field space where the Hamiltonian for perturbations is unbounded from below, due to the cubic interaction (|47p . Unstable 
solutions present ghosts either in the scalar sector (branch II) or in the vector sector (branch I and II), and are represented by 
red. In the branch intersection, there is strong coupling (S.C.) for both the vector and scalar modes, since they do not have 
kinetic terms, but the vector instability remains. 
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IV. MAXIMAL SYMMETRIC SOLUTIONS AND INSTABILITIES FOR ARBITRARY 03 AND 04 



In this section, we extend the previous discussion to the case of arbitrary a^, and a^. Including the contributions 
corresponding to non- vanishing and 04, the Lagrangian in the decouphng hmit reads 



The expression for the X^il is given in eqs p^ . To study this Lagrangian and its maximally symmetric solutions, 
we proceed exactly as in section ITTll performing the shifts of tensor and scalar modes as in eqs ([T3l) - ([T4| . We skip 
the details as they are a straightforward generalization of the arguments of the previous section, and present only the 
final results. 

In order to eliminate the terms that are linear in the fields in the resulting Lagrangian (tadpoles), one obtains the 
following set of two equations that generalize eqs (|24p - (P5)) : 



\q-H^-{1- Co) [2 - Co + Aai{\ - co)^ + (4 - 5co + c^)] = , (58) 
[3- 2co + 12^4(1 -co)2 + 3a3(3- 4co + c2)] =0. (59) 

After diagonalizing the Lagrangian (although at fourth order in perturbations a mixing between h and tt remains, 
when 03 7^ — 4a4), one obtains the following kinetic term for the scalar field 

Ck^n,. = [3-2co + 12a4(l-co)'+3a3(3-4co + c2)]' 

-3i7^ [1 + 3^3(2 - Co) + 12^4(1 - Co)] jTrDTT . (60) 

These results were already obtained in The new piece of information is the Lagrangian for vector perturbations, 
which up to third order in perturbations becomes (obtained following the same method of section fill B I) 



[3 + 12a4(co - If - 2co + 3a3(3 + (co - 4)co)] trFa 

(61) 



4co 

[-1 + 3a3(-2 + Co) + 12a4(-l + Co)] trHtrFj 



4co 

+^ [3 + 9a3 + 12a4 - 4(1 + 6a3 + 12a4)co + 9(^3 + 4a4)co] trHFs . 

The previous Lagrangian contains only the firsts of an infinite series of terms including vectors, but the terms above 
are enough for our purposes. 

Within this system, we would like to determine maximally symmetric configurations and discuss their stability. As 
in the previous section, eq. (|59|) admits two branches of solutions: one corresponding to Minkowski space, and the 
other to a maximally symmetric space-time. These are the generalization of the two branches of solutions discussed 
in section UlI CI Since we are interested in characterising self-accelerating de Sitter configurations, we will not discuss 
the Minkowski branch any further in the main text, relegating its discussion to Appendix [Cj We will, instead, focus 
on the maximally symmetric option 7^ for solving eq. (|59p . The solution for co is then 



^ (1-1- 303) 3 (aa + 4^4) ±\J{1 + ^a^f - 3 (a3 + 4^4) 

'° " 3(a3 + 4a4) ' ^^^^ 

l + 3a3T|a5| 

where, in order to render the formulae less cumbersome, we define (as done in [Toj) 



^ - (l-H3a3)^-3(a3-K4a4) . (64) 



Notice that in the special case of = — 4a4, the positive branch Cq disappears and only Cq^ remains physical. In the 
following, we will impose > in order to have well-defined square root in eq. (j62p . and co > as discussed after 
eq. (j38l) . Upon substituting the previous values for co in the expression for the scalar kinetic terms, we find 

Ckrn,. = ±3|a5|-ff'^n7r. (65) 
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Substituting the value of cq from eq. (l63l) into eq. ([58)) . we find the following expression for 



— Hq + Xq , 



± ^ 1 l + 3a3T2|a5| 
° - 3(l + 3a3T|a5|)2 



(66) 



The results so far suggests that self-accelerating configurations with > may be allowed in this theory, by a 
suitably choice of and a^; this is indeed the result of j9|]. However, let us see what happens to vector fluctuations. 
The coefficient in front of the vector kinetic term can be written in terms of as 



[(1 + 3^3)^ 



4co 



-(co - cd)(co - Cg )tri^2- 



(67) 



The solutions for cq, eq. (j63p . automatically imply that the vector fluctuations are strongly coupled. After substituting 
our solution for cq of eq. (1631) in the Lagrangian (j6ip . we find 



^l«5| 



(l + 3a3T Ictsl 
4(2 + 3a3T|a5 

trntrF2 - 2tr(nF2) 



trntrFa - 2tr(nF2) 



(68) 



This is a generalisation of the result for q;3 = q;4 = 0. This Lagrangian leads to equations of motion that are at 
most second order in time derivatives. However, using arguments identical to the ones developed in section THI C[ it is 
possible to see that, when as ^ 0, the Hamiltonian associated with this system is unbounded from below if one moves 
along certain directions in the moduli space of perturbations. Hence, for the very same reasons explained in section 
If H C[ also these solutions with arbitrary and 04 (and keeping as ^ 0) are unstable when cq satisfies eq. ([63)) . In 
Fig. 2 we summarise the stability of the solutions in the general case. The special case as = is more subtle, since 
both scalar and third order vector Lagrangian vanish, and it is related with the configurations analyzed in [2j]. We 
will postpone its study to the future. 
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FIG. 2: Stability analysis of screening solutions (brancii I) and self-accelerating solutions (branch II) for a general 03 and 
a4, obeying the condition (1 -I- Sas)^ > Q5. Solutions with instabilities are shown in red and healthy ones in blue. The 
precise values of H^, and the points where Ao = depend on the explicit values of 03 and cti. There are two branches of 
self-accelerating solutions (solid horizontal lines), which lack of kinetic terms for vector perturbations (strong coupling), and 
present an instability when cubic interactions are considered. If 03 and as are chosen correctly, one may have expected a 
healthy model of de-Sitter without a cosmological constant Ao, but this is not the case due to the instabilities described in the 
main text. There are three branches of screening solutions with = (solid vertical lines), which are connected at cq — ■ 
In the case of (1 -I- Za^f < a%, the diagram is similar, but the screening solutions are healthy in the intermediate segment 
(between and c^) and unstable elsewhere. Finally, if (1 -I- 803)'^ = a%, or equivalently 03 — —4a4, the two self-accelerating 
branches collapse into a single one and it becomes similar to Fig. [1] 
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V. DISCUSSION 



In this work we studied, in detail, the dynamics of vector fluctuations in the decoupling limit of A3 massive gravity, 
with the main aim of investigating the stability of self-accelerating configurations. The decoupling limit of A3 massive 
gravity propagates 5 physical degrees of freedom around any background: 2 tensors, 2 vectors, and 1 scalar. While it is 
well known that, in the decoupling limit, scalar fluctuations organize so to form the few scalar Galileon combinations, 
the structure of the vector Lagrangian is more subtle. It contains an infinite number of terms, one at each order 
in perturbations. We developed a method that allows to determine, in a fully covariant way, the vector Lagrangian 
at each order in perturbations. Around self-accelerating configurations, the vectors become strongly coupled, being 
characterized by vanishing kinetic terms. This has been shown perturbatively in 9] and non-perturbatively for 
spherically symmetric perturbations in [lOj. In this paper, we showed non-perturbatively that the strong coupling 
of the vector perturbations is an inevitable consequence of selecting the self-accelerating branch of the solutions. 
However, the vector acquire dynamics thanks to contributions at third and higher order in perturbations; at third 
order they form a p-form Galileon combination, which ensures, at this order, that the associated equations of motion 
contain at most two time derivatives. At fourth and higher orders it is not automatically so, and the equations 
of motion contain higher derivatives. However these higher derivatives are harmless, since one can show that field 
redefinitions recast the Lagrangian in a manifestly healthy form, which does not propagate additional degrees of 
freedom. We find interesting to reveal this structure of the decoupling limit Lagrangian when vectors are included. 
It generalises the standard scalar Galileon combinations to p-form Galileons and indicates in, yet another way, how 
the BD 'sixth mode' is not present around special backgrounds. 

However, the strong coupling of the vector degrees of freedom around self-accelerating solutions is problematic, since 
it may turn one of the available five degrees of freedom into a ghost, making these configurations generically unstable. 
A careful Hamiltonian analysis of the system, indeed, reveals that the Hamiltonian is unbounded from below along 
generic directions in the space of fluctuations. It would be very interesting to understand whether such instability 
can be avoided, for example imposing symmetries that prevent to move along dangerous directions in moduli space. 

The strong coupling also means that quantum corrections are not controllable in general, and we should stress 
that our analysis here is purely classical. It is known that the Galileon terms, without a vector, are protected from 
quantum corrections [27j . though we expect that the Maxwell term for vector perturbations will receive quantum 
corrections [o^. It will be important to study quantum corrections when vectors are included. 

In our work we did not only limit our attention to self-accelerating solutions, but also to more general maximally 
symmetric solutions in the presence of a bare cosmological constant. We investigated the two branches of vacuum 
solutions. The first one contains the self-accelerating configurations; it is characterized by a non-vanishing curvature 
, the strong coupling of the vector sector, and aforementioned instability. The second branch is such that curvature 
vanishes, even in the presence of a bare cosmological constant. These solutions realise the screening configurations 
already studied in Q. We confirm that, around these solutions, all fluctuations (vector included) are well behaved, 
provided the bare cosmological constant lies within an interval. 

Our conclusions have been obtained in the decoupling limit, that corresponds to an analysis of the theory in a well 
deflned, but limited, regime of scales. It would be interesting to reproduce and interpret these results from the full 



theory perspective, and point out other instabilities besides the ones analyzed in 11[. Such analysis will be carried 



on in a further publication [25| , aimed also to clarify issues associated with the dynamics of scalar fluctuations in the 
full theory, which were found to exhibit a strong coupling behavior (no kinetic terms) when expanded at quadratic 
order around certain cosmological backgrounds [19|, |20| . 
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Appendix A: Higher derivatives and the scalar-vector Lagrangian 

In this appendix we investigate the strongly coupled scalar vector Lagrangian of eq. (I45p . that we rewrite here: 

= - (trHtrF^ - 2trHF2) - - (trtf - trHtrHF^ + trHFHF) + . . . , (Al) 
6 9 
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where the dots represent higher order perturbations. The vector has no kinetic terms, and its dynamics is controlled 
by perturbations starting at third order (as already noticed in 0, [13, The third order Lagrangian has a 

p-form Galileon structure [l^], and leads to equations of motion that contain at most two time derivatives: 

(F^.^pf- 2 (F^, = 0, (A2) 

F^^"- " n^^ + - F^^^n = o . (A3) 

More delicate and interesting is the dynamics associated with the fourth order Lagrangian . This has not the 
structure of a p-form Galileon (see eq. (1461) ). that in our context would read 

^p-/or™ ^ 1 ^^.^2 _ 2 trHF^) + 6o I ^trF^ trH^ - (trH)^ - 2 trH trHF^ + 2 trH^F^ + trHFRFj , (A4) 

with bo an arbitrary coefficient; then one might wonder whether it leads to additional degrees of freedom. 

We have many independent proofs, obtained following different procedures ^ that show that A3 massive gravity 
does not propagate a sixth ghost mode around any background. In the present context, we can make this fact 
concretely manifest in various ways. For example, by showing that combining the equations of motion, it is possible 
to recast them in such a way that no higher derivatives are present. Or, alternatively, determining a field redefinition 
that, starting from the Lagrangian of eq. (|Aip . allow to recast the 4th order part in a form that coincides with the 
p-form Galileon of eq. (IA4p . so to render the system manifestly healthy. Let us discuss this second approach; the field 
redefinition that we need is 

F F - i n F. (A5) 

By plugging the previous redefinition into the £3 part of (jAip . we completely remove the fourth order part of the 
Lagragian (in other words, we generate the healthy scalar- vector lagrangian as in eq. (|A4p with bo ~ 0). 

On the other hand, we know that our vector-scalar lagrangian contains an infinite number of contributions, while 
until now we only considered perturbations up to fourth order. What happens at higher order can be investigated by 
focussing on a simple, but representative example, in which the total lagrangian can be computed at all orders and 
resummed. We consider a situation in which only the following components of the vector are switched on: 

= {Ao{y), 0, Ayit), 0) . (A6) 

We also suppose, again for simplicity, that the scalar depends only on the (i, x, y) coordinates, and the cross derivatives 
vanish (11;^,^ = if /j, ^ v). Then the action for the vector can be calculated at all orders in perturbations and can be 
resummed exactly (as done in |l3|) finding for cq = 3/2 

f2 n 

Ca = ^^^^ f^j^Y^j 

where Fyt = (dyAo — dtAy). Notice that second derivatives appear at the denominator; such a system would naively 
seem to lead to the propagation of additional modes. This Lagrangian has the same structure of the one discussed in 
for the spherically symmetric case, and more recently in [26] for a purely scalar theory. Since we have only one 
component Fty switched on for the tensor F^^ it is simple in this case to exhibit a field redefinition (that generalizes 
(|A5p at higher orders) that removes the dangerous terms from the Lagrangian. More instructive is in this case to 
analyze the equations of motion, and showing that by manipulating them they can recast in a form that manifestly 
does not lead to the propagation of an additional mode. The equation of motion for Ay tells us that 

(iOyAo-8,Ay)U..A^^ ^ flMi^M^U constant . (A8) 

The equation of motion of tt gives 

„n / (dyAo — dtAy)^ 'nxx\ , . . 1 . /Ar>N 

c/j — 2 — + parts that contam at most two tmie derivs — . (A9) 

\ (3 4- Iltt - 11^) / 

Plugging (jASp into (jA9p . one easily finds that all terms containing three or more time derivatives cancel out. Con- 
sequently, these higher time derivatives terms do not lead to additional 'ghost' mode. This shows directly, at least 
within this special ansatz, that the resummed Lagrangian does not propagate additional degrees of freedom. It would 
be very interesting to extend these results to the complete Lagrangian for the vector-scalar modes, without making 
any assumption on the field profiles. 
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Appendix B: Turning on a background profile for the vector field 

In this section, we turn on a background profile for tlie vector field, and show that perturbations around the 
resulting configuration always contain a ghost degree of freedom. The discussion essentially repeats the analysis of 
[lOll. a lthough in the present set-up we make use of the Lagrangian only up to third order in perturbations (while 
in [10| we considered a resummed Lagrangian in a spherically symmetric set-up). We then focus on the following 
scalar- vector Lagrangian up to third order in perturbations 

j^tMrd ^ _3ij-2^n7r_ 1 [UttF^^F'^'' - 2d^,d„iTF^'PF''p] . (Bl) 

It is not hard to check that the corresponding equations of motion admit the following solution 

sol 2 Qq 2 ITir,\ 

- = , (B2) 

Af = -Qo {x^ + y^ + z^) , (B3) 

this configuration does not break the de Sitter symmetry of the physical metric. Perturbing the third order action 
around this solution, and calling the fields 

Ao = -Qo{x'' +y'^ + z'')+Ao, (B5) 

we obtain (using the divergenceless condition 9tAo — diAi) the following induced kinetic terms at second order in 
the small perturbations tt, Aq 



,2„d order _ 3^2 -q- ^ ^^UA^ + ^Pf^ (B6) 



The last term in the previous expression does not contain time derivatives. Let us now focus on the time derivatives 
contained in the first two terms, which read 



£2nd order ^ 3ij2 ~ _ ^^^^2^^ _^ pj^^^g ^j^j^ derivatives (B7) 



Defining an effective scalar tt = tt — |^ one then finds that the previous Lagrangian contains always a ghost when 
H ^ Q m. agreement with [l^]. While in this appendix we focussed on the case = a^^ = Q it is straightforward to 
extend this analysis to the more general case, obtaining the same conclusion. 



Appendix C: Minkowski branch of the solutions for arbitrary as and 04 



In this appendix, we discuss the solutions to eqs when selecting = to solve eq. ([5^ . In this case, both 

scalar and vector perturbation have non- vanishing kinetic terms, so it is sufficient to study the sign of their coefficients 
to investigate the stability of the configurations. Using defined in (|63p . the kinetic terms can be written as 



27 



(co - c+)(co - Co ) 



ttDtt, 



(CI) 



>Lfcm,A„ = -. (cq - Cq )(co - Co )tri'2, (C2) 

4Co 

where co is a solution of the tadpole equation ()58p . which can be written as 

F(co) EE 2 + 4^3 -f- 40:4 - 3(1 -I- 3a3 -I- 4q!4)co + (1 + 603 + 12a4)cg - (03 + 4a4)cg = Aq. (C3) 

The stability of vector perturbations depend on the solution for cq, which is a cubic equation in co if 0:3 ^ —Aa^. 
The cubic function F{co) has two stationary points at cq = c^, F'(cJ) = where Ao = = —2H^ {H^ is defined 
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in ([66])). Thus there are three branches of sohitions and at the stationary points Cq = Cq, where the two branches of 
solutions meet. These points are also the boundaries of the different stability behaviours. For example, if > — 4a4, 
and as and 04 are such that > 0, then the vector is healthy between cq > Cq and < cq < Cq (see figure [2|). 
On the contrary, if < —Aa^, but > 0, then the vector is healthy only within the region Cq < cq < Cq . For 
as = — 4a4, eq. (|C3I) is satisfied when 



3 + 6a3 ± ^1 + 4 (1 + Sag) Ao 

= ht^) ■ ^^'^ 

The kinetic terms for scalar and vector seem to vanish when = —'ia^; however it is not really the case since 
have (ofa + 4a4) at denominator (see eq. (|62p). Depending on the sign chosen in eq. (|C4p . the kinetic term for the 
vector can then be rewritten as 



r (l + 3a3) v/l + 4(l + 3a3) Ao 

2(3 + 6^3 ± ^1 + 4 (1 + 3^3) Ao) 

while the kinetic term for the scalar has no risk to become negative definite. Focussing for simplicity on the case 
1 + 3a3 > 0, the plus sign in eq (|C4[) leads to a ghost in the vector sector. So only the negative sign is allowed, and 
the allowed solution for cq is 



1 + 2 (1 + 303) - v/1 + 4 (l + 3a3) Ao 



With this choice of sign, we have a well defined square root in (jC6p and a healthy kinetic term in (IC5[) if Ao lies within 
the interval 

^ < Ao < 2 + 3a3 , (C7) 



4 (l + 3a3) 

generalizing the result discussed in Section IIII CI obtained in the case 0(3 = 0(4 = 
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